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Abstract. Cosmic voids are a promising environment to characterize neutrino-induced ef¬ 
fects on the large-scale distribution of matter in the universe. We perform a comprehensive 
numerical study of the statistical properties of voids, identified both in the matter and galaxy 
distributions, in massive and massless neutrino cosmologies. The matter density field is 
obtained by running several independent V-body simulations with cold dark matter and neu¬ 
trino particles, while the galaxy catalogs are modeled by populating the dark matter halos 
in simulations via a halo occupation distribution (HOD) model to reproduce the clustering 
properties observed by the Sloan Digital Sky Survey (SDSS) II Data Release 7. We focus on 
the impact of massive neutrinos on the following void statistical properties: number density, 
ellipticities, two-point statistics, density and velocity profiles. Considering the matter density 
field, we find that voids in massive neutrino cosmologies are less evolved than those in the 
corresponding massless neutrinos case: there is a larger number of small voids and a smaller 
number of large ones, their profiles are less evacuated, and they present a lower wall at the 
edge. Moreover, the degeneracy between ag and is broken when looking at void properties. 
In terms of the galaxy density field, we find that differences among cosmologies are difficult 
to detect because of the small number of galaxy voids in the simulations. Differences are 
instead present when looking at the matter density and velocity profiles around these voids. 
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1 Introduction 

IV-body simulations of structure formation show that the spatial distribution of matter resem¬ 
bles a network: the so-called cosmic web. In this web, knots represent regions of high density, 
known as dark matter halos, connected by filaments, i.e. relatively thin and low density 
regions. Finally, filaments are usually surrounded by large regions with very low densities: 
voids. 

Besides the anisotropies in the cosmic microwave background, the vast majority of the 
constraints on the values of cosmological parameters arise from the observations of only one 
element of the cosmic web: the dark matter halos (or the galaxies residing within them). Much 
attention has been put in understanding the spatial distribution and statistics of galaxies and 
dark matter halos. This is because it is believed that those objects are biased tracers of the 
underlying matter distribution. Thus, by measuring their spatial correlations on large scales 
it is possible to constrain the values of the cosmological parameters. Additional information 
on the value of the cosmological parameters can be obtained by measuring the abundance of 
dark matter halos as a function of their mass, i.e. via the halo mass function (see for example 
[!])• 

Recently, it has been pointed out that cosmic voids can also be used to constrain the 
values of cosmological parameters [2-12], As a different element of the cosmic web, the 
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information embedded into the properties of voids will complement the one obtained from 
the dark matter halos. That additional information can then be used to break degeneracies 
and to further constrain the parameters of the cosmological model. Among the different 
observables that can be used to retrieve cosmological information with voids there are: the 
distribution of void ellipticities [4], the Alcock-Paczynski test [5, 11, 12], the integrated Sachs- 
Wolfe effect [6] and the weak lensing effect [13]. 

However, cold dark matter and baryons are not the only massive particles present in the 
universe that could play a role in the formation of structure, and in particular in the evolution 
of voids: there are also neutrinos. The standard model of particle physics describe neutrinos 
as fermions of spin 1/2, neutral and massless particles. However, neutrino oscillations exper¬ 
iments have pointed out that at least two out of the three neutrino families are massive (see 
[14, 15] for review). These observations are one of the clearest evidences pointing towards the 
existence of physics beyond the standard model. Thus, one of the most important questions 
of modern physics is to know the masses and hierarchy of the three neutrino families. 

The existence of a neutrino cosmic background is one of the predictions of the Big Bang 
theory. The presence of a hot component in the total matter density held affects cosmological 
observables in many different ways. At the linear level massive neutrinos modify the matter- 
radiation equality time and the growth of matter perturbations [16, 17]. Those effects are 
commonly used to put upper limits on the sum of the neutrino masses from a variety of 
cosmological observables such as the anisotropies in the cosmic microwave background (CMB), 
the spatial clustering of galaxies, and the clustering properties of neutral hydrogen absorbers 
in the Lya forest [18-38]. The upper limits on the sum of the neutrino masses obtained from 
cosmological observables m v < 0.12 eV at 95% C.L. [38]) are much tighter than those 
achieved from laboratory experiments. 

Therefore, cosmological observables are currently the best probe to constrain neutrino 
masses and thus, it is timely to investigate the impact of massive neutrinos on the large scale 
structure. The impact of massive neutrinos on different non-linear observables such as the 
matter power spectrum, the abundance and spatial distribution of dark matter halos and the 
galaxy properties has been recently addressed by a number of authors [39-60]. 

There are reasons to expect that massive neutrinos affect more strongly the properties 
of cosmic voids than the properties of dark matter halos. Given their large thermal velocities, 
neutrino clustering within dark matter halos and galaxies will be limited to the neutrinos 
populating the low-momentum tail of the Fermi-Dirac distribution [54], In other words, the 
relative contribution of relic neutrinos to the total mass of dark matter halos is much smaller 
than the cosmic ratio A consequence of this fact is that most of the properties of 

the dark matter halos, such as their mass function or their spatial bias, can be described in 
terms of the CDM density held, rather than the total matter density held [45-48, 57, 58, 61]. 
On the other hand, the large thermal velocities of the neutrinos will prevent their evacuation 
from cosmic voids. This should manifest as an extra mass within voids that will affect their 
overall evolution. It is thus expected that voids in a massive neutrino universe would be 
smaller and denser. This fact was hrstly noted by [7], who studied the properties of voids in 
massive neutrino cosmologies using the Lya forest. 

In this paper we investigate for the hrst time the impact of massive neutrinos on cosmic 
voids, as identihed in the spatial distribution of galaxies and in the underlying matter density 
held. We study the properties of voids in the spatial distribution of matter by running large 
box A r -body simulations in cosmologies with massless and massive neutrinos. The spatial 
distribution of galaxies is modeled using a halo occupation distribution (HOD) model and 
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requiring that the resulting mock galaxy catalogues reproduce, for a given galaxy population, 
the observed number density and two-point correlation function of a particular survey. We 
use the publicly-available code VIDE [62] to identify the voids in both the matter and galaxy 
distribution. The void catalogs obtained with this void finder - and watershed void finder in 
general, to be discussed below - depend on the resolution of the N-body simulation considered. 
However, in this work we are interested in the relative differences between cosmologies with 
and without massive neutrinos, rather than absolute results. Therefore, we compare catalogs 
and void properties obtained from simulations at the same resolution. 

This paper is organized as follows. In Sec. 2 we describe the N-body simulations we have 
run for this project. The method used to identify voids in both the matter and the galaxy 
fields is briefly outlined in Sec. 3. A comparison between the different properties of voids 
detected in the matter density held in the different cosmological models is shown in Sec. 4. 
In Sec. 5 we present the properties of voids identified in the galaxy distribution. Finally, we 
draw the main conclusions in Sec. 6. 

2 N-body simulations 

We ran the V-body simulations using the TreePM code GADGET-III [63]. The simulation 
parameters are summarized in table 1. The values of the cosmological parameters, for all 
simulations with the exception of the L0s8 and of the L0s8CDM runs, are the ones found 
by the Planck collaboration [64]: Q m = 0.3175, tit, = 0.049, = 0.6825, h = 0.6711, 

n s = 0.9624, A s = 2.13 x 10 -9 . We have run simulations for four different cosmological 
models: a model with massless neutrinos and three models with massive neutrinos (three 
degenerate species) with masses ^2m v = 0.15 eV, = 0-3 eV and 'Yh rn v = 0.6 eV. In 

massive neutrinos cosmologies we set H c d m = H m — Hb — £l v where Q u h 2 = ^ m v / (94.1 eV). 
We chose these values for the neutrino masses because we want to study and isolate the effect 
of massive neutrinos on voids properties. In simulations with the sum of the neutrino masses 
lower than current tightest bounds < 0.12 eV [38] the effect of massive neutrinos 

may not have been seen properly or it would have been largely contaminated by cosmic 
variance. Note further that some groups are claiming that concordance cosmology between 
Planck CMB data and cluster abundance or weak lensing could be achieved by allowing a 
non-zero neutrino mass of 0.3 eV (±0.1 eV) (see however the discussion in the latest Planck 
cosmological parameter paper [64] ). In order to investigate the m u — as degeneracy, we 
have also run simulations with massless neutrinos, but with a value of as equal to the one of 
the ^2 m u = 0.6 eV model. The as value can be computed from the total matter or from the 
CDM power spectrum; we denote the simulations with the former as L0s8, the ones with the 
latter as L0s8CDM. 

In this paper we investigate the properties of voids as identified in the matter and 
galaxy density fields, in cosmologies with massless and massive neutrinos. Our simulations 
can thus be divided in two different categories: 1) simulations with low resolution and large 
volumes, used to detect voids in the matter distribution and 2) simulations with relatively 
high resolution and low volumes, used to isolate voids in the spatial distribution of galaxies. 

The low resolution simulations consist of ten independent realizations, for the same 
cosmological model, each of them following the evolution of 256 3 CDM particles (and 256 3 
neutrinos for massive neutrinos cosmologies) within a periodic box of size 1000 Mpc /h. We 
chose to work with such a low number of particles because it is difficult to run the void finder 
VIDE on more particles, due to the large memory required. The high resolution simulations 
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Name 

Box Size 
(Mpc/h) 

(eV) 

V c dm 

N v 

e 

(kpc/h) 

C8 

(z=0) 

realizations 

L0 

1000 

0.00 

256 3 

0 

100 

0.834 

10 

L0s8 

1000 

0.00 

256 3 

0 

100 

0.693 

10 

L0s8CDM 

1000 

0.00 

256 3 

0 

100 

0.717 

10 

L15 

1000 

0.15 

256 3 

256 3 

100 

0.801 

10 

L30 

1000 

0.30 

256 3 

256 3 

100 

0.764 

10 

L60 

1000 

0.60 

256 3 

256 3 

100 

0.693 

10 

HO 

500 

0.00 

512 3 

0 

25 

0.834 

1 

H15 

500 

0.15 

512 3 

512 3 

25 

0.801 

1 

H30 

500 

0.30 

512 3 

512 3 

25 

0.764 

1 

H60 

500 

0.60 

512 3 

512 3 

25 

0.693 

1 


Table 1. Specifications of our N-body simulation suite. The first letter of the simulation name 
indicates whether it is a low-resolution (L) or high-resolution (H) simulation. The values of the 
following cosmological parameters are the same for all the simulations (with the exception of the 
simulation L0s8 in which A a = 1.473 x 10~ 9 and L0s8CDM where A s = 1.576 x 10 -9 ): = 0.3175, 

= 0.049, fl A = 0.6825, h = 0.6711, A s = 2.13 x 10" 9 , n s = 0.9624. 


only contain one single realization, which traces the evolution of 512 3 CDM particles (and 
512 3 neutrinos) in a box of 500 Mpc /h size. In all simulations the gravitational softening is 
set to 1/40 of the mean linear inter-particle distance. 

We generate the initial conditions of the simulations at z = 99 by displacing the particle 
positions, that are initially located in a regular grid, making use of the Zel’dovich approxi¬ 
mation. In simulations with massive neutrinos we use the massive neutrinos transfer function 
and a mass weighted average between the CDM and baryons transfer functions to set up the 
initial conditions of the neutrino and the CDM particles, respectively. In massless neutrinos 
simulations we use the total matter transfer function for setting up the initial conditions of 
the CDM particles. Notice that, in massive neutrino cosmologies, the displacements and the 
peculiar velocities of the CDM and the neutrino particles in the initial conditions have been 
generated taken into account the scale-dependent growth present in these cosmologies. 

3 Void Finder 

We identified voids using the publicly available Void Identification and Examination (VIDE) 
toolkit [62], which uses a substantially modified version of ZOBOV [65] to perform a Voronoi 
tessellation of the particles and then a watershed transform to group the Voronoi cells into 
a hierarchical tree of subvoids and voids [66]. A final catalog of voids is then built using the 
criteria that the voids must be larger than the mean particle separation, so for our analysis we 
will only consider voids larger than 4.0 Mpc jh in the CDM field (low resolution simulations) 
and voids larger than 10 Mpc/h in the galaxy distribution (high resolution simulations). 
While voids naturally form a nested hierarchy, we only consider top-level (i.e., parent) voids 
in this work. 

When identifying voids in the matter density field, we run VIDE on top of the cold 
dark matter particles only, even if we are considering cosmologies with massive neutrinos. In 
principle we should select voids in the total matter field, including also neutrinos. However, 
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VIDE cannot discriminate between two different particle populations having two different 
masses. Therefore we must select voids in one of the two density fields and the CDM is the 
one mainly responsible for the evolution of the cosmic structures. It would be interesting to 
understand the differences arising in defining voids in the CDM or in the total matter field 
and which of the two is the fundamental field to identify voids. When selecting voids in the 
galaxy distribution we instead run VIDE on top of the whole mock galaxy catalogue. 

The output of VIDE is a catalog which contains much information about the identified 
voids. In this work we use only the following: position of the center, particle members, 
effective radius (i? e ff) and elhpticity. The effective radius R e g is defined as the radius of a 
sphere containing the same volume as the watershed region that delimits the void, and the 
center is the volume-weighted center of all the Voronoi cells in each void. 

4 Void in the matter distribution 

In this section we present the analysis for voids identified in the matter density field, which is 
modeled via our low resolution iV-body simulations. As explained in the above section, voids 
are selected in the matter density field by running VIDE on top of all CDM particles, for 
every realization of each cosmological model. Our goal is to investigate the impact of massive 
neutrinos on the formation and evolution of cosmic voids. Thus, here we study some of the 
main properties of voids and examine how these depend on the cosmological model. The 
void properties analyzed are: number densities, ellipticities, two-point correlation functions, 
density profiles, and velocity profiles. 

4.1 Number density 

We compute the number density of voids at z = 0 and z = 1 as a function of their effective 
radius R e g, and show the results in Fig. 1. Different colors indicate different cosmologies 
while error bars represent the scatter around the mean value obtained from the 10 indepen¬ 
dent realizations divided by -v/IO. The bottom panel of each plot displays the ratio between 
the number density of voids in massive to massless neutrinos cosmologies. We find that in 
cosmologies with massive neutrinos the abundance of small voids is larger than in massless 
neutrinos cosmologies, whereas the number density of big voids is highly suppressed in mas¬ 
sive neutrino cosmologies. These trends take place both at z = 0 and z = 1, while relative 
differences are smaller at z = 1. 

These results are in perfect agreement with our expectations: neutrinos have large ther¬ 
mal velocities, which gives rise to two main effects. On one hand it avoids their clustering 
within dark matter halos, and on the other hand it makes them less sensitive to void dynam¬ 
ics. Therefore, as a first approximation, the neutrinos density contrast can be approximated 
as 5 U = 0 and the matter density as <5 m = 5 c dmPcdm/Pm + SuPv/Pm — ^cdmPcdm/Pm, where the 
subscripts ’cdm’ and’m’ stand for CDM and matter, respectively. This brings an extra mass 
inside the voids that will slow down their overall evolution. It is thus expected that voids in 
a massive neutrino universe would be smaller and denser, and therefore appear less-evolved. 
Moreover, the higher the neutrino mass the higher the ratio p u /p m and the additional mass 
inside the voids. Therefore the relative differences with respect to ACDM increases with the 
neutrino mass. We also expect differences to become smaller at higher redshift. At higher red- 
shifts voids will be denser overall, and the additional mass given by neutrinos will play a less 
important role and the void number density will be more similar among different cosmologies. 
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Figure 1 . Number density of voids as function of their effective radius. The top panels show 
the results for different cosmologies: red line is the ACDM cosmology, purple, blue and green lines 
correspond to Yl rn v = 0-15 eV, Y^ 171 " = 0-3 eV and Y2 m v = 0-6 eV cosmologies, respectively. The 
bottom panels show the ratio between the void number density in massive and massless neutrino 
cosmologies. Left and right plots display results at redshift z = 0 and 2 = 1, respectively. 


We conclude that the number density of voids, identified in the matter density field, is 
very sensitive to the neutrino masses. 

4.2 Ellipicity 

We now investigate the distribution of voids ellipticities in cosmologies with massive and 
massless neutrinos. For each cosmological model we compute the fraction of voids with a 
given ellipticity and show the results in Fig. 2 at redshifts 2 = 0 (left) and 2=1 (right). 
The error bars represent the scatter around the mean value obtained from the 10 independent 
realizations normalized by \/l0. The bottom panels show the ratio between the models with 
massive neutrinos to the model with massless neutrinos. We find differences of the order of a 
few percent between ACDM and 0.15 eV cosmologies for void with ellipticities e < 0.4, where 
statistical error bars are relatively small. The differences increase as the sum of neutrino 
masses increases and it reaches 20% for the 0.6 eV cosmology for both redshifts 2 = 0 and 
2 = 1. We find that relative differences among models slightly decrease at 2 = 1. 

With our interpretation that voids in massive neutrino cosmologies appear younger than 
those in massless neutrino cosmologies, it is straightforward to understand the results. Voids 
in cosmologies with massive neutrinos will be in a earlier evolutionary stage and therefore we 
would expect to find more voids with low ellipticities than in a ACDM cosmology. For the 
same reason, there is a deficit of voids with large ellipticites in massive neutrino cosmologies. 

4.3 Correlation function 

The structure of the void correlation function has been well studied (for example, in [67]), 
which identified the void-exclusion scale and the effects of galaxy bias. Here we investigate 
the clustering properties of voids in the different cosmological models. For each cosmology, 
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Figure 2. Number of voids with a given ellipticity, normalized by the total number of voids. The 
top panels show results from different cosmologies: red line is the ACDM cosmology, purple, blue 
and green lines correspond to = 0.15 eV, Y2 Tn v = 0.3 eV and ^ m„ = 0-6 eV cosmologies. 

The bottom panels show the ratio between the same quantity in a massive and a massless neutrino 
cosmologies. Left and right plots display results at redshift z = 0 and z = 1, respectively. 




R [h 1 Mpc] R [h 1 Mpc] 

Figure 3. Two-point correlation function of voids with an effective radius in a given range. Solid 
lines show the results for the ACDM cosmology whereas the dashed lines indicate the results for the 
Y2 rriv = 0.6 eV cosmology; different colors are related to different bins in the void radii as reported 
in the legend. Left and right plots display results at redshift z = 0 and 2 = 1, respectively. 


we split our void catalogues in different groups depending on their size R e g and compute the 
correlation function of the void centers belonging to the same group. The correlation function 
is measured using the Landy-Szalay estimator [68]. The random catalogue contains a number 
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of points 20 times larger than the void catalogue. 

Figure 3 shows the results at redshift z = 0 (on the left) and z = 1 (on the right), 
with error bars given by the scatter around the mean value obtained from the 10 independent 
realizations divided by vTO- For clarity we show the correlation functions for only two 
cosmologies: the solid lines refer to ACDM and the dashed ones refer to ^ m v = 0.6 eV 
cosmology. Different colors indicate different ranges in R e g. The results show that, for both 
cosmologies, the two-point correlation function goes to zero on large scales. It increases as 
R decreases and it reaches a maximum, whose height and position depend on the size of the 
voids considered. Finally, going towards smaller distances, the correlation function decreases 
until it reaches the lowest boundary -1. 

The behavior of the correlation function on small scales is due to exclusion effects: top- 
level voids are extended objects that cannot overlap. In fact, if we imagine them as spheres 
of radius R e g, the probability of finding two void centers at a distance smaller than the sum 
of their effective radii (exclusion scale) is zero, which corresponds to a value of the correlation 
function equal to -1. We would naively expect to have a sharp transition from the positive 
correlation to the exclusion regime. However, this is not the case for several reasons. First 
of all, the presence of different void sizes in each of the considered groups. Secondly, the 
void finder used in this paper does not return spherical voids, but under-dense regions with 
a complicated shape, and therefore it can happen that the distance between two void centers 
is smaller than the sum of their two effective radii. It can be seen that the exclusion scale 
increases with the effective radius of the voids considered. 

The presence of a positive peak in the two-point correlation function arises from two 
different processes: the rise of the non-linear clustering of voids going to smaller scales and 
the exclusion effect. We can also understand why the clustering of voids in massive neutrino 
cosmologies is lower with respect to the massless neutrino case if we take into account that 
voids in the former are younger than those in the latter. Indeed, if we compare the clustering 
of voids at z = 0 and at z = 1, we find that the clustering decreases with redshift. Thus, since 
voids in the massive neutrino cosmologies are effectively younger, they are expected to exhibit 
a lower amplitude in the 2pt-correlation function, as our results show. At higher redshifts 
the properties of voids in cosmologies with massive and massless neutrinos get closer, and 
therefore it is natural that differences in their clustering properties decrease. 

4.4 Density profile 

Here we investigate the impact of massive neutrinos on the shape of the void density profiles. 
We compute the cold dark matter, massive neutrino, and total matter density profile in voids 
of different sizes. 

The profiles have been computed in the following way. For each simulation box, we select 
the voids whose effective radius he within a certain interval. For each void, we compute the 
density held inside shells around the center, then we stack all the shells of the different voids 
to obtain the average density profile at a given distance from the void center. We repeat 
the procedure for all the ten realizations with the same cosmology, then we average over 
the realizations. The error that we present is the scatter between the different realizations 
normalized by \/l0. We repeat this procedure for different intervals in the void radius and 
for all the cosmologies considered in this paper. 

Figure 4 shows the average cold dark matter (solid lines) and neutrino (dotted lines) 
density profiles around void centers with radius R e s- 10-11 Mpc /h (top), -R e ff=16-18 Mpc/h 
(center), and R e g= 20-25 Mpc /h (bottom). Left and right plots show results at redshift z = 0 


- 8 - 






R [h~ l Mpe] JJI/T 1 Mpc] 




R [h -1 Mpc] 
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Figure 4. Average cold dark matter (solid lines) and neutrino (dotted lines) density profiles around 
voids with different sizes: i? o g=10-ll Mpc /h (top), J? e ff=16-18 Mpc/h (center), i? e ff=20-25 Mpc /h 
(bottom). Left and right panels show results at redshifts z = 0 and z = 1, respectively. Red, purple, 
blue and green lines show the 0.0, 0.15, 0.3 and 0.6 eV cosmologies, respectively. The dashed black 
lines indicate the mean value of the void radii in the selected range and two times the same quantity. 


and z = 1, respectively, and different colors indicate different cosmologies. Let us focus first 
on the cold dark matter density profile. We can immediately notice that small voids present 
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a compensated wall at the edge, whereas the large ones have a negative density profile, as 
already shown in [69]. Moreover, fixing the cosmology, voids present a typical behavior: 
large voids are emptier; i.e., they present deeper underdensities in the core than small ones. 
However, we have checked that this behavior can change depending on the resolution of the 
considered simulations. 

For voids with a given radius, we find that the height of the wall becomes higher when 
decreasing the neutrino mass. On the other hand, the inner part of the density profile is 
emptier in the massless neutrino case and it becomes denser as the neutrino mass increases. 
Voids evolve in time building a higher wall and becoming progressively emptier and emptier. 
Again we can understand the difference between the cosmologies: the impact of the neutrino 
mass can be described as slowing down the evolution of the void profiles and giving a less 
evolved universe in these regions. 

The neutrino profile is found to follow the corresponding cold dark matter one. Around 
very small voids, the neutrino profile presents an overdensity at the void wall, whereas it 
presents troughs around very large voids. Any departure from the mean background density 
(both overdensities and troughs) is greater for larger neutrino masses. This corresponds to 
the neutrinos with lower thermal velocities, which are more susceptible to the presence of the 
cold dark matter gravitational held. 

In Fig. 5 we show the total matter density profile, which is the most relevant quantity 
since most of the observables depend on the total matter distribution, e.g. weak lensing and 
the ISW effect. We present the results for the same void sizes as above and at z = 0 (left) and 
z = 1 (right). The bottom panel of each plots shows the ratio between the density profiles of 
the three massive neutrino cosmologies with respect to ACDM. The main features observed 
in the cold dark matter profiles are overall preserved, but enhanced here. In the void core, the 
differences between ACDM and the 0.15 eV cosmology are at the level of 1-3 %, depending 
on the void radius, and they reach 10% for the 0.6 eV case. Near the wall and for the very 
small voids, the difference is at the 5 % level for 0.15 eV and at the 20% level for 0.6 eV 
cosmologies. All the differences are more pronounced at z = 1. 

4.5 Velocity profile 

We now present the radial velocity profiles of the cold dark matter, neutrino, and total matter 
fields. The first two have been computed directly from simulations in the following way. We 
select voids with radii in a given range and then we compute the radial velocity profiles of 
particles inside a shell of radius r around each void center, using 

1 N r- 

< r ) = > 
i— 1 

where N is the number of particles in the shell and f) are the coordinates of the particles 
with respect to the void center. Next, we stack all the shells from the different voids for 
computing the average velocity at a given distance from the void center. Finally, we average 
over the ten realizations with the same cosmology. The error associated with the profiles is 
the scatter between the ten realizations divided by \/l0- The velocity of the total matter field 
is instead computed as the density weighted average between the velocity profiles of the cold 
dark matter and neutrino particles, via 

Vcdm{r)pcdm(r) + Vy(r)p u {r) 

Pcdm(r) + p u (r) 


(4.2) 


Vm (r) 


(4.1) 
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Figure 5. Average total matter density profiles around voids with different sizes: i? e ff=10-ll 
Mpc//i (top), i? e ff=16-18 Mpc /h (center), and i? e ff=20-25 Mpc /h (bottom). Left and right panels 
show results at redshifts 2 = 0 and 2 = 1, respectively. Red, purple, blue and green lines show the 
0.0, 0.15, 0.3 and 0.6 eV cosmologies, respectively. At the bottom of each panel we display the ratio 
between the results from the massive neutrino cosmologies and the ACDM one. The vertical dashed 
black lines indicate the mean value of the void radii in the selected range and two times the same 
quantity. 
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where the subscripts m, cdm and v stand for matter, cold dark matter and neutrinos, respec¬ 
tively. The associated errors are computed via error propagation. Other velocity estimators 
have been used in literature, e.g. the Voronoi weighting estimator [69]. However, we have 
checked that the relative differences between cosmologies (massive vs massless) are also re¬ 
produced by this estimator. 

The results for the cold dark matter (solid lines) and the neutrino (dotted lines) fields 
are shown in figure 6. As for the density profile, we consider three ranges for the void radius 
and we plot on the left/right the results at redshifts z = Q/z = 1. 

The cold dark matter particles have a positive radial velocity inside the effective radius 
(first vertical dashed lines from the left), meaning that the inner part of the voids is expanding 
and becoming more and more empty. The behavior outside the effective radius changes 
depending on the void size. Around very small voids (top panel), the cold dark matter velocity 
becomes negative, meaning that the particles are moving towards the void. Therefore, there 
is an ongoing construction of the wall around the edge of the void. The void will eventually 
shrink and later collapse because it is surrounded by an overdense region (this is the so-called 
void-in-cloud effect). Around bigger voids, the radial velocity is always positive; it presents a 
peak near half of the void radius and it decreases on large distances from the center. Therefore, 
the outer region of the void expands slower than the inner part, producing nevertheless a 
concentration of mass around the edge of the void. Looking at different cosmologies, the 
average radial velocity is higher in ACDM and it decreases as the sum of the neutrino masses 
increases. Going from redshift z = 0 to z = 1, the velocities increase for all void radius and 
cosmologies and differences between models become more pronounced. 

As expected, we find that the neutrino radial velocity field follows the cold dark matter 
one. The velocity is positive inside the void radius and remains positive also outside, apart 
from the case of very small voids, where it becomes negative. However, the neutrino velocity 
appears to be smaller than cold dark matter one, which is only due to cancellation effects. 
Neutrinos have large thermal velocities that make them free-stream in every direction, there¬ 
fore the average velocity is expected to be close to zero. To better understand this effect, 
Fig. 8 shows the average measured positive and negative neutrino radial velocity profile, which 
correspond to the profiles for the outgoing and incoming neutrinos, and their theoretical pre¬ 
diction. The predicted mean thermal velocities of incoming and outgoing neutrinos are equal 
due to spherical symmetry and it is computed as Vi n = v ou t — 160/2(3eV/Sm„) km/s [16]. 
The two measured velocities are similar for small voids and small neutrino masses, but they 
differ otherwise. It is the gravitational interaction with the void and the surrounding matter 
which gives rise to this and to a positive (or negative) average neutrino velocity. At higher 
redshift (z = 1), the neutrinos are faster and they feel less the dynamics of the underlying 
cold dark matter structure. Indeed, the average neutrino radial velocity profile is lower at 
z = 1 than at z = 0. 

The results for the total matter radial velocity field are presented in Fig. 7. The main 
features described for the cold dark matter radial velocity also apply here; however, the 
differences between different cosmologies are more pronounced. In the bottom panel of each 
plot we show the ratio between the results for the massive neutrino cosmologies to the ACDM 
case. At redshift z = 0 and in the inner part of the voids there are differences at the order 
of 5-10-15% for 0.15-0.3-0.6 eV cosmologies, respectively. Larger differences can be seen in 
the outer parts of the voids, in the ranges i? e g=10-ll Mpc /h (top) and i? e ff=16-18 Mpc jh 
(center). At redshift z = 1 the differences are slightly more pronounced. 
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Figure 6. Cold dark matter (solid lines) and neutrinos (points) velocity profile measured from 
simulations around voids with different sizes: i? o ff=10-ll Mpc /h (top), i? e ff=16-18 Mpc /h (center), 
and i? e ff=20-25 Mpc /h (bottom). Left and right panels are computed at redshift 2 = 0 and 2 = 1, 
respectively. Red lines indicate the ACDM cosmology, and purple, blue and green lines show the 0.15, 
0.3 and 0.6 eV cosmologies. The dashed black lines indicate the mean value of the void radii in the 
selected range and two times the same quantity. 


4.6 fijy-og degeneracy 

We also investigate the well known degeneracy between Q, v and erg [46, 56, 59]. For instance, 
the effects of massive neutrinos in many observables such as cluster number counts, can be 
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Figure 7. Total matter velocity profile measured from simulations around voids with different sizes: 
i? e ff=10-ll Mpc /h (top), 16-18 Mpc /h (center) and 20-25 Mpc /h (bottom). Left and right panels 
display results at redshift z=0 and z=l, respectively. Red lines indicate the ACDM cosmology, purple, 
blue and green lines show the 0.15, 0.3 and 0.6 eV cosmologies. At the bottom of each panel there are 
the residuals between the massive neutrino cosmologies and the ACDM one. The dashed black lines 
indicate the mean value of the void radii in the selected range and two times the same quantity. 


14 - 

























































R [IT 1 Mpc] 

Figure 8. Radial velocity profile of outgoing (solid lines) and incoming (dashed lines) neutrinos 
around voids with different sizes: i? o g=10-ll Mpc//i (top-left), 16-18 Mpc /h (top-left), and 20-25 
Mpc /h (bottom). The solid horizontal lines show the theoretical mean radial velocity; purple, blue 
and green indicate the 0.15, 0.3 and 0.6 eV cosmologies. At the bottom of each panel there are the 
ratios between the velocity of outgoing and the incoming neutrino for each cosmology. The dashed 
black lines indicate the mean value of the void radii in the selected range and two times the same 
quantity. 


mimicked by a cosmology with massless neutrinos but with a lower value of <7g, which is 
the r.m.s. value of the linear fluctuation in the mass distributions at 8Mpc//i. We want to 
understand if this is the case for voids too. 

We use the simulations in ACDM with erg taken from the m u = 0.6 eV cosmology 
(L0s8 and L0s8CDM) to identify voids in the matter density field. If the degeneracy is present 
also in voids, these simulations should give catalogs very similar to the one of the 0.6 eV case. 
Moreover, we consider the two different erg values because it has been shown that for halos 
the degeneracy arises when using the CDM field, since it is the fundamental one in driving 
the evolution of overdense structures [46]. We must check if this is the case for underdense 
structures as well. 

In order to test all of this, we first compare the number density of voids and we plot the 
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Figure 9. Number density of voids as function of their effective radius. The top panel shows 
the results for different cosmologies: red and green lines are the ACDM cosmology with <t s equal to 
the one in the ^ m „ = 0.6 eV case computed from the total matter and the CDM power spectra, 
respectively. The blue line correspond to the ^ m v = 0.6 eV cosmology. The bottom panels show 
the ratio between the void number density in massless and massive neutrino cases considered. 


results in figure 9. The ACDM cosmology with erg computed from the total matter field (red 
line) matches quite well the predictions from the 0.6 eV case up to void radii equal to ~17 
Mpc /h. Above that, the departure is significant, giving a lack of big voids, and it reaches 
a disagreement at the level of 70-80% for voids with size ~30 Mpc /h. Things are different 
for the ACDM cosmology with as computed from the CDM field (green line). In this case 
there are fewer small voids and this translates into a worse match with the massive neutrino 
cosmology, which reach the 5% for very small voids. On the other hand, the number density 
of big voids increases, mitigating the discrepancies at large R e s, which is here at the level of 
~30%. 

In both cases, the results for the m v = 0.6 eV cosmology cannot be reproduced, which 
means that the Q u -as degeneracy is not present in the number density of voids. Voids, unlike 
halos, do not present this degeneracy because they are sensitive to different regions of the 
linear power spectrum. To be more precise, big voids should be influenced by scales in the 
P(k) much larger than the ones important for the halos. On these large scales, the amplitude 
in the ^ m u = 0.6 eV case is higher than in the two considered massless neutrino cosmologies, 
and more power corresponds to a higher number density of voids. 

We also perform the same comparison using the density and velocity profiles. The results 
are plotted in figure 13, which we placed in Appendix A. 

5 Voids in the galaxy distribution 

It has to be noticed that matter is not (unfortunately) directly observable and usually we rely 
on large scale structure tracers like galaxies. The distortion in the shape of galaxies caused 
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by the weak gravitational effect can be used to constrain cosmology. The weak lensing signal 
will depend on the total matter density distribution and it would be interesting to answer 
the following two questions: 1) what is the matter density profile of the voids identified in 
the galaxy held? 2) what is the extent to which massive neutrinos impact the reconstructed 
profile? Recently some work have been done to address the first questions (see for example 
[70-72]); here we focus on the second one. Therefore, in this section we present the analysis 
for voids identified in the galaxy distribution, which is obtained by post-processing IV-body 
simulations. First, we explain how we populate halos with galaxies in order to generate a 
mock galaxy catalogue. Then we show the main results regarding the study of galaxy voids. 
Since the volume of our high-resolution simulations is relatively small and the number density 
of galaxy-voids is low, we cannot perform the analysis as for the voids in the cold dark matter 
held, studying all the void properties in different cosmologies. Therefore, we present only the 
number density of galaxy voids and the study of matter inside the identified galaxy voids. In 
particular, we show the matter density and velocity profiles around these objects. 


5.1 Halo Occupation Distribution 

In order to construct a mock galaxy catalog, we need to populate with galaxies the dark 
matter halos in IV-body simulations. A commonly used tool is the so-called Halo Occupation 
Distribution (HOD) model. In this framework, the distribution of galaxies with respect 
to dark matter halos is described by a few parameters that can be calibrated in order to 
reproduce some particular features of the observed galaxy population. Here we calibrate the 
HOD parameters to reproduce the number density and the two-point clustering statistics of 
galaxies in the main sample of the Sloan Digital Sky Survey (SDSS) II Data Release 7 [73]. 

The HOD model requires two ingredients: 1) the probability distribution p(N\M) of 
having N galaxies inside a halo of mass M and 2) the way in which galaxies positions and 
velocities are related to the underlying matter particles. In our HOD model the first ingredient 
is described by three parameters: M m j n , a and M\. We assume that halos with masses below 
Mj ni n do not host any galaxies, whereas halos with masses above M m - m host one central galaxy 
and a number of satellite galaxies following a Poissonian distribution with a mean equal to 
(M/Mi)". We use the following equations: 


(N C \M) 


1 if M > M min 
0 if M < Mmin 


(N S \M) 


(M/Mi)" 

0 


if M > M m i 
if M < M m j 


(5.1) 


Our second HOD ingredient states that the central galaxy sits in the center of the correspond¬ 
ing halo and that the distribution and velocity of the satellites follow exactly the ones of the 
underlying cold dark matter particles inside the halo. This means that the galaxy bias and 
the velocity bias with respect to the cold dark matter held are both equal to 1. The values 
of the HOD parameters obtained for galaxies with magnitudes M r — 51og 10 h = —21.0 in the 
four different cosmologies are shown in table 2. 

We focus our study at redshift z = 0, where there are good measurements of the clus¬ 
tering of different galaxy populations. Now that we have created a mock galaxy catalog for 
each cosmology, we can run the void-finder VIDE on top of the galaxy held and identify the 
galaxy voids. 


5.2 Number density 

For each cosmology we compute the number of galaxy voids as a function of their size R e g and 
we show the results in Fig. 10. On the upper panel we plot the cumulative number of galaxy 
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E 777 -*' 

(eV) 

Mi 

(M e /h) 

a 

Afmin 

(M e /h) 

0.0 

1.22 x 10 14 

1.38 

4.92 x 10 12 

0.15 

1.17 x 10 14 

1.38 

4.79 x 10 12 

0.3 

1.18 x 10 14 

1.48 

4.50 x 10 12 

0.6 

9.86 x 10 13 

1.46 

4.16 x 10 12 


Table 2. Values of the HOD parameters, for four different cosmologies and for galaxies with magni¬ 
tudes M r — 51og 10 h = —21.0. 



Figure 10. Cumulative fraction of galaxy voids at z = 0. The top panel shows results from different 
cosmologies: red line is the ACDM cosmology, and purple, blue and green lines correspond to 0.15, 0.3 
and 0.6 eV cosmologies. The bottom panel shows the ratio between the same quantity in a massive 
and a massless neutrino cosmologies. 


voids normalized by the total number and the associated errors, which are computed using 
error propagation and assuming a Poisson distribution in R e g- for the galaxy voids. Different 
colors indicate different cosmologies. On the bottom panel we display the ratio between the 
results for the massive neutrino cosmologies and the ones for the ACDM model. We can 
notice some differences among the cosmologies only for voids with R e g > 30 Mpc /h, where 
anyway the errors are very large. Moreover, it is interesting to notice that the size of the 
voids identified in the galaxy field are larger than the ones found in the cold dark matter 
field. This is just a consequence of the number density of the tracers used to identify voids: 
a smaller value for the number density of tracers brings to larger voids. 

5.3 Density profile 

Here we compute the total (cold dark matter + neutrinos) matter density profile in the 
identified galaxy voids, once they are divided into four groups, depending on their effective 
radius: i? e ff=20-25, 25-30, 30-35, 40-45 Mpc jh. Figure 11 shows the profiles (top panel of each 
plot) only for 0.0 eV and 0.6 eV cosmologies, while the residuals between massive and massless 
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Figure 11. Density profile of the total matter field (cdm + neutrinos) around galaxy voids with 
different sizes: R e g = 20 — 25 Mpc /h (top-left), R e g = 25 — 30 Mpc jh (top-right), R e g = 30 — 35 
Mpc /h (bottom-left), and R e g =40 — 45 Mpc /h (bottom-right). In each plot, the main panel shows 
the result for the 0.0 eV (red) and the 0.6 eV (green) cosmologies. The bottom panels display the 
ratios between the results from the massive neutrino cosmologies (0.6 eV but also 0.15 eV in purple 
and 0.3 eV in blue) to the ACDM one. The vertical dashed black lines indicate the mean value of the 
void radii in the selected range and two times the same quantity. 


cases (bottom panel of each plot) for all the massive neutrino cosmologies are displayed. As 
in the previous section, the vertical dashed lines indicate the effective radius and two times 
the same quantity. 

We immediately notice that the matter profile here, in galaxy voids, is very different 
from the one around voids identified in the cold dark matter field (see figure 5). First of all, 
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galaxy voids are less empty. Secondly, there is a compensated wall even around very large 
voids. These differences are due to the fact that galaxy voids incorporate within them many 
small cold dark matter voids, together with small halos and filaments. Their edge presents 
an overdensity in the galaxy field, which can sit only in very big overdensities in the cold 
dark matter field. In other words, this should be due to the galaxy bias b [74]. If we assume 
that voids are patches where the enclosed number density contrast of the tracers is equal to 
a certain threshold 5t = —0.8, and that the tracers are galaxies with a density field S g = b5 m 
biased with respect to the matter field, then the corresponding enclosed matter density will 
be S m = St/b. If the galaxy bias is b ~ 2, the enclosed matter density will be 2 times denser 
then the galaxy one. 

Focusing on the comparison among cosmologies, we can notice the same trend as in 
the cold dark matter voids: the matter profiles of ACDM are emptier and they present 
a higher wall. However, here the small galaxy void statistics makes this trend less clear. 
Moreover, the differences among ACDM and massive neutrinos cosmologies are more noisy, 
but not less pronounced. In small voids there are smaller differences, whereas in large voids 
the differences are of the order of 10-20% in the inner part of voids and for all the massive 
neutrino cosmologies; the discrepancies are less pronounced near the compensated wall. 

5.4 Velocity profile 

We present here the analysis for the velocity profile of matter around galaxy voids. We select 
the galaxy voids in the same size ranges as for the density profile. We compute the radial 
velocity of matter as a density weighted average between the velocity profiles of cold dark 
matter and neutrino fields, as shown in equation (4.2). The associated errors are computed 
via error propagation, using the dispersion around the mean values of density and velocity 
profiles of cold dark matter and neutrinos obtained from the stack of the selected voids in the 
simulation. 

The results are shown in Fig. 12. Again we plot the profiles for ACDM and 0.6 eV 
cosmologies only, but we present the residuals with respect to ACDM for all the massive 
neutrino cosmologies considered in this paper. The main feature to notice is the particular 
shape of the profiles. They are positive in the inner part of the galaxy voids, they reach zero 
around the effective radius and they are negative outside. Even the profiles around very big 
voids share the same shape, contrary to what happens in cold dark matter voids (see figure 7). 
It seems that in these regions the cold dark matter is building a wall near the effective radius, 
with the additional effect of shrinking the underdensity. 

If we compare the cosmologies we notice that inside galaxy voids the velocity is higher 
in ACDM, whereas in outer regions the trend is not unique for all the void sizes. However, for 
most of the cases, the massive neutrino cosmologies have slower motion towards the center of 
the galaxy voids. 

6 Conclusions 

The aim of this work is to analyze the relative differences in the void properties due to the 
presence of massive neutrinos. We have run two sets of IV-body simulations, one with low 
resolution and another with higher resolution. Each set of simulations have been computed 
in four different cosmologies, i.e. with neutrino masses equal to 0.0 eV, 0.15 eV, 0.3 eV and 
0.6 eV. We have used the low resolution simulations to identify voids in the CDM particle 
distribution and the high resolution ones for studying voids in the galaxy field. 
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Figure 12. Radial velocity profile of the total matter field (cdm + neutrinos) around galaxy voids 
with different sizes: R e g = 20 — 25 Mpc /h (top-left), R e g = 25 — 30 Mpc /h (top-right), R e g = 30 — 35 
Mpc /h (bottom-left), and R e s = 40 — 45 Mpc /h (bottom-right). In each plot, the main panel shows 
the results for the 0.0 eV (red) and the 0.6 eV (green) cosmologies. The bottom panels display the 
ratios between the results from the massive neutrino cosmologies (0.6 eV but also 0.15 eV in purple 
and 0.3 eV in blue) to the ACDM one. The vertical dashed black lines indicate the mean value of 
the void radii in the selected range and two times the same quantity, whereas the colored numbers 
indicate the number of stack voids in each simulation. 


For the voids in the CDM held, we have initially studied the number density of voids in 
the four different cosmologies at redshift z = 0 and z = 1. In both cases, the number density of 
small voids is higher in massive neutrino cosmologies, whereas the number density of big voids 
is lower. The difference between ACDM and cosmologies with massive neutrinos increases 
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with the neutrino mass, but it decreases with the redshift. This can be understood in terms 
of the total mass contained in the void. Indeed, neutrinos have high thermal velocity which 
prevents them to feel the void dynamics and to create a deep underdensity around the void 
center. This translates into an extra mass within the void that will evolve slowly. Therefore 
voids in the massive neutrino universe are smaller. However, the additional mass has a smaller 
effect at high redshift, where voids are denser, thus resulting in smaller differences at z = 1. 

Secondly we studied the distribution of void ellipticities. In ACDM we found that the 
fraction of voids with low ellipticity is smaller than in the massive neutrino cases. On the 
other hand, the fraction with high ellipticity is larger. This can be interpreted in term of 
evolutionary stages of voids in the different cosmologies: voids in massive neutrinos universes 
are younger than in their massless neutrino cosmology corresponding model. 

We computed the correlation function of voids having radius within specific ranges. The 
results showed the exclusion effect on very small scales, which is due to the fact that voids are 
extended objects that cannot overlap. This effect, in addition to the presence of void clustering 
on small-intermediate scales, creates a positive peak, which is placed at different distances 
depending on the void size considered. The height of the peak changes with the cosmology: 
it increases with lower neutrino masses. Looking at the differences between redshift z = 0 
and z = 1, we can understand that this is due to the fact that voids in massive neutrino 
cosmologies are effectively younger. Moreover, differences between cosmologies become less 
important with redshift. 

Another interesting property of voids is their density profile. We computed the CDM, 
the neutrino, and the total matter density profiles for different void sizes and in the four 
different cosmologies. We found that small voids present a compensated wall at the edge, 
whereas large voids have a non-positive profile. Again, we can understand the differences 
between massive and massless cosmologies once we consider that voids are effectively less 
evolved in the presence of massive neutrinos. Indeed, voids evolve by evacuating particles, 
becoming progressively emptier and building the wall. For each void size in both CDM and 
matter density field, ACDM cosmology presents higher walls and emptier cores, in agreement 
with this explanation. Instead, the neutrinos profile is flatter, given the neutrinos’ high 
thermal velocities. They also tend to follow the CDM one: around small voids there is 
a positive overdensity of neutrinos in correspondence with the wall, whereas the neutrino 
profile presents usually an underdensity in the core of big voids. All departures from the 
mean background get more pronounced as the neutrino mass increases. Focusing on the total 
matter profile, which is the most important one since most of the observables depend on the 
total matter distribution, we obtained the following results. In the void core, the differences 
between ACDM and the 0.15 eV cosmology are at the level of 1-3%, depending on the void 
radius, and they reach the 10% for the 0.6 eV case. Near the wall and for the very small 
voids, the difference is at the 5% level for 0.15 eV and at the 20% level for 0.6 eV cosmologies. 
The differences increase at redshift z = 1. 

The last investigated property of voids identified in the CDM field is the radial velocity 
profile of CDM, neutrino, and total matter around voids. The CDM and total matter profiles 
have the following behavior: the radial velocity is positive inside the void radius around which 
it reaches the maximum, then it decreases and it becomes negative for very small voids. It 
remains positive, but at a small value, for large ones. This means that particles inside voids 
move towards the edge, evacuating and expanding the voids, whereas the outer particles have 
different behaviors depending on the void size. In small voids they move towards the center 
participating in the building up of the wall and eventually making the void collapse (the 
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void-in-cloud effect). Instead, in big voids they go far away, but at somewhat smaller velocity 
than the ones in the inner parts, giving rise to a concentration of mass around the edge of the 
void. The neutrino radial velocity follows the same behavior, but it has a smaller magnitude, 
and this is due only to a cancellation effect. Indeed, neutrinos have high thermal velocities 
that make them free-stream in every direction and the average velocity is close to zero. 

We also studied the degeneracy between cr 8 and f \ v running two set of 10 simulations in 
ACDM with o '8 taken from the m u = 0.6 eV cosmology. We saw that this degeneracy is not 
present when focusing on void properties since their number density, but also their density 
profile and velocity, is different in massless and massive cosmologies with same er 8 . 

In the second part of this work we have analyzed the high resolution N-body simulations. 
We have populated these simulations with galaxies, following the HOD prescriptions, and we 
have identified voids in the galaxy distribution. We have shown the cumulative fraction of 
voids at redshift z = 0. There are some differences among the cosmologies only for voids 
with R e ff > 30 Mpc /h, where the errors are very large due to the small volume probed by 
our high resolution simulations. Moreover, the voids are bigger than the ones selected in the 
CDM field, since the number of tracers here is low given the small box-size. Current surveys 
can span a larger volume: for example the Sloan Digital Sky Survey 1 III (SDSS-III) with the 
Data Release 9 (DR9) CMASS sample can covers an effective volume of nearly 1.5 (Gpc//r) 3 , 
which is about 20 times bigger than our box-size. 

Finally we studied the total matter density and velocity profiles around the galaxy voids, 
since future surveys like Euclid 2 and DESI are expected to significantly increase the number 
of observed galaxies. Furthermore, the weak lensing signal will allow to measure the total 
matter density field. The matter density profiles present some differences with respect to the 
ones in the CDM voids: they are less underdense and they present the wall even around very 
big voids. We argued that this is probably due to the galaxy bias with respect to the matter 
density field. Instead, for what regards comparison among cosmologies, we found the same 
behavior as in the CDM voids. The velocity profiles show a peculiar shape shared by all the 
void sizes analyzed, despite of what happens for the CDM voids. The profiles are positive in 
the inner part of the galaxy voids, they reach zero around the effective radius and they are 
negative outside. In the inner regions the velocity is higher in ACDM, whereas in the outer 
part of the voids the trend is not unique. 
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A Gj,-ct 8 degeneracy: profiles 

We present the second part of the analysis for investigating the £l u -ag degeneracy using voids. 
Here, we focus on the comparison between the ^ m„ = 0.6 eV cosmology and the massless 

1 http://www.sdss.org 

2 www.euclid-ec.org 

3 desi.lbl.gov 
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Figure 13. Left panel: matter density profile. Right panel: radial velocity profile. In both cases 
colored lines show the profiles for the Y^ rn v = 0.6 eV cosmology, whereas black lines indicate ACDM 
cosmology with the same of the massive neutrino case, when it is computed using the total matter 
density field. The bottom panels show the ratios between the results from the two cosmologies. 


case with same <7s computed from the total matter linear power spectrum. Figure 13 shows 
the density and velocity profiles of the matter held around voids selected in ranges of R e g. 
The colored lines show the results for the massive neutrino case and black lines indicate the 
massless one. The mismatch in the density profile is around 2-4% in the inner part of voids 
and it increases in the outer part, whereas for the velocity profiles it reaches the 10% in the 
void core. 

B Void-by-void comparison among different cosmologies 

In section 4.4 we have compared the density profiles of voids with the same fixed size. It is 
also interesting, mainly for model building purposes, to study the differences among different 
cosmologies by performing a void-by-void comparison. The procedure we use to perform a 
void-by-void comparison is as follows. We first select all voids of a given size in the ACDM 
cosmology. Then, for each of the selected void, we find the corresponding one in a massive 
neutrino cosmology by searching the void whose center he closest to the center of the ACDM 
void. Finally, we use these voids to compute the mean density profile in the massive neutrino 
cosmology. We have explicitly checked that this simple procedure works well when dealing 
with large voids, whereas for small voids the comparison is not always satisfactory. Therefore, 
we limit our analysis to big voids. In figure 14 we compare the density profiles for CDM, 
neutrinos and total matter for voids in ACDM and )T) m v = 0.6 eV cosmologies. For the 
latter cosmology we use both procedures, i.e. voids with radii within a given range (green 
lines) and voids selected by matching them with those of the ACDM cosmology (blue line). 

The comparison between the density profiles of voids in the ACDM cosmology and the 
corresponding voids in the massive neutrino cosmology confirms that the presence of massive 
neutrinos makes the voids matter evacuation slower and therefore they are less empty in the 
center and they present a lower wall at the edge. The comparison between void profiles in 
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the same cosmology m v = 0.6 eV) but found using the two different ways (blue and green 
lines) is enlightening. Indeed, voids selected using the ACDM catalogue appear to be smaller 
than the one selected looking at the void size, since their slope is closer to the center. This 
is even more clear when plotting the void radii in the void-by-void selection for ACDM and 
s f2,m v = 0.6 eV cosmology (lower row in figure 14). The black line indicates the equality 
between the radii in the two cosmologies and most of the points are below that line, which 
means that majority of voids in massive neutrino cosmology are smaller that the ones in the 
massless case. This is even stronger evidence of what we expected: the presence of massive 
neutrinos slow down the evolution (expansion) of voids. 
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Figure 14. Comparison void-by-void between the ACDM and the ^ m„ = 0.6 eV cosmologies. 
We have selected all voids in the ACDM cosmology with radii i? e fj=16-18 Mpc /h (left column) and 
l?eff=20-25 Mpc /h (right column) and we plot their density profiles in red. The upper panels show 
the CDM and neutrino density profiles whereas middle panels display the results for the total matter 
density profiles. For the 0.6 eV cosmology we have chosen the voids using two different procedures: we 
select all voids in above radii range (green lines) and we choose the voids by matching their positions 
to these of the ACDM cosmology (blue lines). The vertical dashed black lines indicate the mean value 
of the void radii in the selected range and two times the same quantity. The lower panels show the 
void-by-void comparison between the radii of the ACDM and the 0.6 eV voids. Each pixel in those 
panels displays the number of void pairs as a function of the void radius i?o.OeV in the ACDM (x-axis) 
and the radius of the corresponding void i?o.6eV in the 0-6 eV cosmology (y-axis). The solid black 
lines represent the curve i?o.6eV = -Ro.OeV- 
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